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The stationary cosmological model without closed timelike curves of Godel 
type is obtained for the ideal dust matter source within the framework of the 
teleparallel gravity. For a specific choice of the teleparallel gravity parameters, 



o 

O 

this solution reproduces the causality violating stationary Godel solution in 

general relativity, in accordance with the teleparallel equivalent of general rel- 
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ativity. The relation between the axial-vector torsion and the cosmic vorticity 



is clarified. 
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I. INTRODUCTION 

Cosmological models with rotation attract considerable attention since the early papers 
of Godel [1]. However, there are two fundamental difficulties of the original Godel rotating 
model which made its physical significance problematic: the absence of expansion and the 
presence of the closed timelike curves. The latter is equivalent to the possibility of violation 
of causality in such a spacetime. Much of the research on this subject [2-16] was aimed at 
removing both drawbacks and thereby in constructing a realistic cosmology with rotation 
and expansion. An overview of the corresponding developments and a comprehensive list of 
the relevant references is given in [17]. 

An alternative approach to gravitation is the so called teleparallel gravity (see [18-20] 
where the basic references on the relevant literature can be found) which corresponds to 
a gauge theory for the translation group based on the Weitzenbock geometry [21]. In this 
theory, gravitation is attributed to torsion [22], which plays the role of a force [23], whereas 
the curvature tensor vanishes identically. Such a spacetime (Weitzenbock manifold) is char- 
acterized by a nontrivial tetrad field as the fundamental structure, and this field gives rise to 
the metric as a by-product. The translational gauge potentials appear as a nontrivial part 
of the tetrad field, and they induces on spacetime a teleparallel structure which is directly 
related to the presence of the gravitational field. 

In this work we address the above mentioned problems of a rotating world in the frame- 
work of the teleparallel gravity theory. The general teleparallelism model is characterized 
by the three coupling constants, and as a result the structure of the corresponding field 
equations is less rigid than that of the general relativity theory. Using this fact, one can 
expect that there exist classes of teleparallel models which are free of one or both difficulties. 

More specifically, in the current paper we demonstrate that the completely causal cos- 
mological Godel type solutions do exist in the appropriate teleparallel gravity models. 

The structure of the paper is as follows. In Section II, we review the main features of 
general teleparallel gravity. In Section III, we obtain the teleparallel Godel type solution. 



Conclusions are presented in Section IV. 

II. TELEPARALLEL GRAVITY 

In the teleparallel gravity, the spacetime is represented by the Weitzenbock manifold W4 
of distant parallelism. This gravitational theory naturally arises within the gauge approach 
based on the group of the spacetime translations. In simplified terms (see [24-27] for the 
advanced discussion of gauge theory of gravity in the framework of nonlinear connections), 
this approach introduces at each point of the spacetime manifold a tetrad (coframe) field 
h a ^. Its inverse (frame) h a u satisfies the condition h a M h a u = 8 V . The tetrad field induces 
a teleparallel structure on the spacetime which is directly related to the gravitational field, 
and the Riemannian metric arises as 

9„u = Vab h% h\, (2.1) 

with the Minkowski matrix 77^ = diag(+l, —1,-1, —1). A parallel transport of a tetrad h a M 
between two neighbouring points is encoded in the covariant derivative 

V u h% = d u h% - T% u h\, (2.2) 

where T a ^ v is the Weitzenbock connection. Imposing the geometric condition that the 
tetrads are parallelly transported in the Weitzenbock space-time, we obtain 

V h a = d h a -T a h a = 

This is the condition of the absolute parallelism, or teleparallelism. It is equivalent to 

rv = h a a d u h\ (2.3) 

which gives the explicit form of the Weitzenbock connection in terms of the tetrad. The 
antisymmetric part of connection, 

rpp = yp _ fp (o 4) 



is the torsion tensor of the Weitzenbock connection. The curvature of the Weitzenbock 
connection vanishes identically as a consequence of absolute parallelism. 

It has been shown in [22] that the most general action in the teleparallel gravity coupled 
to the matter has the form 
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— (L + A) + L mat 
Ak 



(2.5) 



S — Id x h 

with h = det o a M , and 

L = Cl T Xtv T Xtv + c 2 T Xtu T utX + c 3 T°\T\ a . (2.6) 

Here c±, C2, and c 3 are the three dimensionless coupling constants, A is a cosmo logical term, 

and k is the Einstein gravitational constant. For the specific choice of the parameters, 

11 . . 

ci = 4, c 2 = -, c 3 = - 1, (2.7) 

the teleparallel gravity reduces to the so called teleparallel equivalent of the general relativity. 
The Lagrangian (2.6) may be rewritten as 

L = S x ™T Xtu , (2.8) 

where 

S =c{T + — [T -T ) + — {g T a - g T a ) . (2.9) 

In the present paper, we use the traditional tensor formalism. Those readers, who prefer 
a more advanced exterior calculus approach, should address [28] for a general discussion of 
teleparallel gravity in the metric-affine approach. Here, we only mention that the gravity field 
Lagrangian can be rewritten in terms of the quadratic contractions of the three irreducible 
torsion pieces with the corresponding coupling constants 01,02,03. Then one can verify 
that the original coupling constants Ci, 02,03 are related to the new parameters as follows: 
c\ = (2oi + a 3 )/3, C2 = 2(02 — Oi)/3, and C3 = 2(03 — Oi)/3, see [28]. 

By performing variation of the action (2.5) with respect to h a fJi , we get the teleparallel 
gravitational field equation 



^d a (hSr) + S U ^T\ X -±5 p x L = kT/ + A5 p x , 



(2.10) 



where T\ p = h a \ 5(hL milt )/(h5h a p ) is the energy-momentum tensor of matter. In the rest of 
the paper we confine our attention to the case of the ideal fluid with 



T x p =(e + p)u x u p -p8 p x . 



f2.11) 



Here e is the energy density, p is pressure, and -u M is the fluid 4-velocity (normalized by 



u^u^ = 1). 



III. THE TELEPARALLEL GODEL TYPE SOLUTION 



The line element of the Godel type cosmological model is given by the interval 



ds 2 = dt 2 - 2 x /^a(t)e mx dtdy - a(t) 2 (dx 2 + k e 2mx dy 2 + dz 2 ), 



(3.1) 



where m, a, k are the constant parameters and a(t) is the time-dependent cosmological scale 
factor. A general analysis [5,17] of the kinematic properties of such a spacetime shows that 
the closed timelike curves are absent in this manifold when k > (a > by definition). 
Using the relation (2.1) we obtain the tetrad: 

(I -^ae mx 0\ 
a 

a\fa^ke mx 



and its inverse is 



h a , 



\0 

/ 1 








h M 








k+c 







a-vV+fc 



a J 

\ 






V a- 1 J 

Using the Lorentz transformation for the tetrad h a ^ 



(3.2) 



(3.3) 



h a n — A. a b h M , 



(3.4) 



we obtain the Weitzenbock connection 



rv = r V + h a (A-^Afc h%. 



(3.5) 



We assume that the Lorentz matrix depends on the cosmological time only. Hence (denoting 
the time derivative by the dot) d v A a c = b l v A a c and consequently 



(A- 1 d v A) b c = SlL b c , 



(3.6) 



with an antisymmetric matrix (Lbc = — L c b) of the form 

( vi v 2 v 3 \ 

Vi -S 3 S 2 



(3.7) 



V 2 S 3 -Si 

\v 3 ~«2 Si / 

From Eqs. (3.2), (3.3) and (3.7), we can now construct the Weitzenbock connection. Its 
nonvanishing components read: 



, a t ( la 

L t t = v 2 \l — — , T xt = a\vi + s 3 



r\ t = v 2 



a 



1 tt — — , J- xt — L zt — — i J- yt — ~ e 



V 0-+ fc 

ae mz A; 

Vo' + fc 
a 

a 



a + /W ' 



, r* 2t = a[v 3 -s 



a 



\ v * + *7 ' 

(s 3 Va + fc + wi x/o-) , 



(3. 



r* 2t = -v\ t = s 2 , r»„ = V2 



ia/ct + fc 



r^ t = s 3 



i , , r% t = — v 2 w — — , r y zt = -si , 

y/a + k a \a + k Va + k 



? ! 3 



Here, as before, the dot denotes a derivative with respect to the time t. 

Assuming that the matter is co-moving with the space, i.e., taking the velocity u a = 5", 
we find the components of the energy-momentum tensor (2.11): 



T t l = e, T x x =T y y = T z z = -p, T y t = - y /Bae mx {e+p). 



(3.9) 



A. Solving the field equations 

In this paper we will study the case of the stationary world with the constant scale 
factor a = const. In this way we will be able to perform a direct comparison of the solutions 
obtained with the Godel solution in the general relativity theory. 

The value of the parameter 2c\ — c 2 is essential. When this is equal zero, i.e., 2c\ = c 2 
[note that this case includes the teleparallel equivalent of general relativity (2.7)], then the 
difference of the ( y v ) and ( z z ) components of the field equation (2.10) yields: 

In other words, we have three possibilities: either m — 0, then the cosmic vorticity vanishes, 
or the coupling constant c 2 = 0, but then also C\ — and the resulting Lagrangian is 
unphysical. Finally, when both m and c 2 are nontrivial, we should have k = — a/2. This 
is exactly the value of the parameter k of the original Godel model which contains the 
closed timelike curves. Accordingly, we will not consider this case as it basically leads to the 
qualitatively same model. 

Thus we instead assume that 1c\ — c 2 7^ 0. Then the ( 2 ')-th component of the field 
equation yields s 2 = 0, whereas ( x z )-th component yields 



s.=vj^- k . (3.11) 



Furthermore, subtracting ( y y ) from ( z z ), we find 



/ 0- . m 2kd + (k + a)c 2 
\l a + k a (2ci - c 2 ) v / oyo" + A: 

Now, the component (/) minus ( x x ) plus the product of {e~ mx / a^/a) times (</) yields the 
equation 



ka{2 Cl + c 2 + c 3 )[k(vf + vj) + (k + a)vl] + v 1 ^Q l + 2^ g. 



2a(k + a) 2 
Here we have denoted the following combinations of the constants: 
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- = 0. (3.13) 



ka(2c 1 + c 2 + c 3 )[2(k-a)c l + 3{k + a)c 2 ] 
2ci - c 2 
P x = [ aCl - (k + a)c 2 ][A(2k 2 + 2ka + a 2 ) c\ - 2o{2k + a) c ± c 2 

-2kac 1 c 3 -2{k + (r) 2 cl-a(k + (r)c2C 3 \/{2c 1 - c 2 ) 2 . (3.15) 

The ( t x )-th and the (t y )-th components of the field equation read, respectively: 

(2 Cl + c 2 )[2(fc + tr)ci - (3fc + 2a)c 2 ] - [2fcci + (k + <t)c 2 ]c 3 = 
U2m 2 v ^ v / ^T^(2c 1 - c 2 )a 2 ' [ ' } 

- a(2 Cl + c 2 + ca)[fc(^ + v 2 ) + (k + a)v 2 } - v^Q 2 + 2^ P 2 = 

2v^e ma; a(A; + CT) 2 ' l ' j 

Here we have introduced further abbreviations: 

g 2 = |(2ci + c 2 ) [2fe(fe + 3(t)ci - (3fc - cr){k + a)c 2 - (k + <x) 2 c 3 " 

+ [4A; Cl + 3(fc + a)c 2 ]c 3 a\/(2 Cl - c 2 ), (3.18) 



p _ (2ci + c 2 + c 3 )[2q-ci + (k + a)c 2 ][(k + a)c 2 - crci] , , 

2 ~ (2ci - c 2 ) 2 • [ ' j 

Let us analyse the two equations (3.13) and (3.17) which can be considered as the alge- 
braic system for v\ and 1*3. Then discarding the denominators, from the difference and the 
sum of the numerators in (3.13) and (3.17) we find 

= _ 2m Pi + kP 2 = 2m [a a - (k + a)c 2 ] 
a Qi - kQ 2 ak(2ci- c 2 ) 

a(2ci + c 2 + c 3 )[k(v 1 + v 3 ) + (k + a)v 2 \ 



a? Qi - kQ 2 
2m 2 P [crci - (fc + a)c 2 ] 
a 2 k(2 Cl -c 2 ) 2 



(3.21) 



Here we denoted 



P = Ak(2k + 3a) c 2 + 2<t(2£; + a) c x c 2 - a 2 c x c 3 

-(k + a) {2k -a)c 2 2 + 2o{k + o) c 2 c 3 . (3.22) 

Note that although the expressions (3.14), (3.15), and (3.18), (3.19) are quite nontrivial, the 
final result (3.20)-(3.21) is rather simple. Substituting (3.20) into (3.21), we obtain explicitly 



9 k + a 

vi 



-3 



A; 



2 2m 2 (2A; + a) [{k + a)c 2 - <JCi] (2 Cl + c 2 ) 
■u 2 H 



(3.23) 



a 2 ka(2ci + c 2 + c 3 )(2ci - c 2 ) 

Obviously, we can put v 2 = from the equation (3.16). Otherwise (3.16) determines the 
parameter k in terms of a and ci, c 2 , C3. 

Finally, taking into account (3.11), (3.12), (3.20), and (3.23), we can write down the 
diagonal components of the field equation (2.10): 

3m 2 A;(2ci - c 2 ) 



2aa 2 
m 2 k(2ci — c 2 ) 



ke + A, (3.24) 

-Kp + A. (3.25) 



2aa 2 

In order to compare with the old Godel model, we now specialise to the case of the dust 
cosmological matter, p = 0. Then equation (3.25) determines the parameter k/a in terms 
of the coupling and geometrical constants, whereas (3.24) yields the energy density of the 
stationarily rotating world: 

m 2 k(co — 2ci) 

e = ^-= — -. 3.26 

naa 2 

Evidently, for the positivity of energy, the teleparallel coupling constants should satisfy 
the inequality c 2 > 2c\. The important point is that the ratio k/a must also be then 
positive. Hence the resulting rotating universe does not contain closed timelike curves, and 
consequently the causality is not violated. 

B. Finding the Lorentz matrix factor 

In order to complete the description of the exact solution obtained, it remains to find 
the Lorentz matrix A a b which enters the tetrad (3.4). It should be noted that from the 
point of view of the resulting Riemannian geometry, this Lorentz factor is irrelevant since 
the spacetime interval (3.1) is invariant under the local Lorentz transformations of the 
tetrad. However, since we deal with the Weitzenbock geometry, the Lorentz factor is the 
essential element of the definition of the tetrad (3.4) and, accordingly, of the connection 
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(3.8). Assuming that the elements of the unknown matrix A a b depend only on cosmological 
time t, we can rewrite the differential equation (3.6) as 



A\. 



K a r L c 



b ■ 



(3.27) 



L\ 



(3.28) 



Fixing the index a, we thus find that the four elements of each row A° a , or A 1 a , or A 2 a , or 
A 3 a satisfy the system of the linear differential equations (3.27) with the constant coefficients 

/ Ui v 3 \ 
v l -s 3 
s 3 -Si 

\v 3 Si / 

The elements of this matrix were determined above in the equations (3.11), (3.12), (3.20) 
and (3.23) and they depend on the coupling constants of the model and on the geometric 
parameters of the ansatz. However, the explicit form of these relations is not needed for 
our purposes. The solution of such a system is straightforward: We assume the standard 
ansatz of the A b a = e xt a b a . Then the constant coefficients a b a are nontrivial, provided the 
characteristic equation is satisfied: 



0. 



-A vi v 3 

v\ —A s 3 

—S3 —A si 

v 3 — si —A 

Evaluating the determinant, we find the algebraic equation 



(3.29) 



A 4 + A 2 (s 2 + s 



v, - v^ 



{s x vi + s 3 v 3 y 



0. 



(3.30) 



It is easy to see that of the four roots X(A), A = 1, ... ,4, two are real and two are purely 
imaginary. When A is a solution of the equation (3.30), the components «% and a a i represent 
the independent free parameters, whereas the rest of the components are expressed as 



a 2 



^ - Vl ~ Vl n a X 2 +S 2 3 ~V 2 i 

(X o Ot ^ = 

ViS 3 -V 3 S x " SiS 3 + ViV 3 



a 1 . 



(3.31) 
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The general solution of the system (3.27) has the form 

A a 6=E C '(A)«%eV)* ? (3.32) 

i=A 

where the C/a) are arbitrary constants. It is thus interesting to note that whereas the 
geometry of the solution is stationary, with the Riemannian curvature and the Weitzenbock 
torsion tensors being independent of the cosmological time, the teleparallel tetrad has a 
nontrivial time dependence via the Lorentz matrix (3.32), in general. 

IV. DISCUSSION AND CONCLUSION 

In this paper we have obtained the teleparallel version of the stationary rotating cos- 
mological solution of Godel type within the framework of the 3-parameter model with the 
general teleparallel Lagrangian (2.6). In particular, we have demonstrated that even for the 
ideal dust matter source there exists a solution without closed timelike curves and thus the 
causality is not violated in such a spacetime. 

For a specific choice of the teleparallel gravity parameters, this solution reproduces the 
causality violating stationary Godel solution in general relativity, in accordance with the 
teleparallel equivalent of general relativity. 

Previously it was conjectured that the axial torsion is directly related to the spacetime 
rotation, showing that such an assumption is valid for the case of the Kerr geometry [29]. It is 
interesting to check whether this is also true for the rotating cosmology under consideration. 

Using the decomposition of the torsion into three irreducible parts under the group of 
global Lorentz transformations [22] , we can straightforwardly calculate the axial- vector part 

A » = * a-p-t (4.1) 

o 

Here rj^ vpa is the completely antisymmetric Levi-Civita tensor. Explicitly, we find from (3.8): 



A '-°- Am =-h(«-Hkh)' A "-fa7FT=^' (4 - 2) 



2 / \ a \ m \ a , 
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For the exact solution obtained, using the equations (3.11) and (3.12), we prove that the 
first three components vanish A 1 = A x = A y = 0, whereas 

AZ = _ m(2c 1 + c 2 )(a + 2k) 
3a 2 (2ci — c 2 )y / o : v / ^ + o" 

As we can see, the axial vector torsion is, indeed, related to the spacetime vorticity. In 
particular, the axial torsion is apparently directed along the axis of the cosmic rotation and 
its value is proportional to the magnitude of the vorticity. However, this relation should not 
be overestimated. Let us recall that k = — a/2 in the original Godel model. Then we find 
from (4.4) that the axial torsion vanishes completely, N 1 = 0, for this case. The existence 
of the rotating cosmological models with the vanishing axial torsion shows that one cannot 
reduce the description of the cosmic vorticity to the properties of the irreducible parts of 
the Weitzenbock torsion tensor. 

The earlier studies of the teleparallel gravity have revealed the one-parameter family of 
models which are compatible with the observational data [22,19,31,32]. Another character- 
istic feature of that family is the existence of black hole solutions which seem to be absent 
in the generic case [28]. On the other hand, exactly this family appears to be plagued with 
the problem of ghost and negative-energy modes as well as of the "hidden" gauge transfor- 
mations [19,30-32], see also a recent study in [33]. In this sense, it is interesting that the 
nontrivial spatial axial- vector torsion (4.1) of our solution prevents the occurrence of the 
special "hidden" gauge transformation [31] of the torsion. At the same time, it is worth- 
while to note that the one-parameter family admits only the causality violating cosmological 
solutions with rotation, whereas the timelike closed curves are absent in the generic case. 
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